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8 SUMMARY 


The condition for plastic buckling of a long flat plate under 
combined shear and longitudinal compression is obtained by using the 
theory of plastic buckling developed for simple loading, provided 
the .ratio of shear to ltsngitudinal compression is assumed constant 
during the loading process. The plate may be elastically restrained 
along the edges. Combinations of shear and longitudinal compression 
were computed for three simply supported plates of 2kS- J Tb aluminum 
alloy, reduced to stress— ratio form, and compared with the known 
interaction curve in the elastic region. Departures from the elastic 
interaction curve were found to be slight, provided the coordinates of 
the interaction curve were modified to allow for variations in moduli. 

INTRODUCTION 

- The combinations of shear and direct stress in the elastic range 
that will cause buckling in an infinitely long flat plate were pre- 
sented in reference 1. The combinations proved to be given with high 
accuracy by the formula 

E 0 + E s 2 = 1 U) 

where R c is the ratio of compressive stress when buckling occurs 

under the combined loading to the compressive stress when buckling 
occurs under compression alone and R a is the ratio of shear stress 

when buckling occurs under the combined loading to the shear stress 
when buckling occurs under shear alone . The interaction curve 
(a parabola) represented by equation (1) was shown to be largely 
Independent of any rotational elastic restraint along the edges of • 
the plate. 
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In the plastic region, solutions have already "been obtained for 
the buckling of an infinitely long flat plate in compression alone 
(reference 2) and in shear alone (reference 3)* The purpose of the 
present paper is to obtain the corresponding solution for any combi- 
nation of compression and shear when applied together. Since the 
■differential equation for plastic buckling in this case has nearly 
the same form as the corresponding equation for the elastic buckling 
of an orthotropic plate , it might be expected that the solutions would 
be similar in the two cases. Balabukh, whose work is discussed in 
reference 4, concluded that f or_ elastic buckling of an ortho tropic 
plate equation (l) still applies with good accuracy. The solution to 
the present problem was sought, therefore, in the same general form 
as equation (l). 

Plasticity theory in its present state is unable to deal with 
problems in combined stress in a general way. However, Ilyushin has 
' shown (reference 5) that if combined stresses are applied in such a 
way that their ratio is a constant during the loading process, the 
combination of stresses may be treated as though it were a -single 
stress. Thus, the theory of plastic buckling already derived is 
applicable to the present problem provided the ratio of shear to 
compression is constant during loading. This condition is often met 
or approximated in practice. The use of the deformation theory of 
plasticity in buckling problems has been justified by Batdorf 
(reference 6). 


SYMBOLS 


a x applied direct stress, positive in compression 

t applied shear stress 


E elastic modulus 

E a secant modulus at stress intensity cr^ 

E^ tangent modulus at stress intensity 

coefficient In buckling formula, a function only of 
length-width ratio of plate and elastic restraint 
in the elastic range and also of stress in the plastic 
range 


stress intensity 


fax 2 + 3t 2 


k 
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b 

h 

D 


plate width 
plate thickness 

tending stiffness of plate in elastic range 


Eh~ 

> 9 


D‘ = 


Eg hr 


S 

e 

X 

0 

Subscripts 

c 

s 

pc 

Pa 


a i 


stiffness of restraining medium 

magnitude of elastic restraint along parallel edges 


/ 4Sb 

W - 


half wave length of buckles 

angle of nodes with normal to plate edges (see fig. 1) 

compression 

shear 

pure compression 
pure shear 

under combined stress 


INTERACTION CURVE 


In the plastic range, the critical stress for pure compression 
is shown to be of the form given in the equation preceding equa- 
tion (28) of reference 2: 


(a x ) pc k pc 


* 2 D f 


= k. 


b c h 


pc 


^ a ^pc rt 2 P 


E 


b 2 h 
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Similarly, the critical stress for pure shear is shown to he of the 
form given in equation (12) of reference 3: 


^ Es ^ps Tt 2 I> 


ps 2 
h n 


In these expressions the subscripts pc and ps refer to pure 
compression and pure shear, respectively. The values of k depend 
upon the aspect ratio of the plate., the conditions of edge restraint, 
and the stress when the plate is stressed into the plastic region. 

When compressive stress a x and shear stress t are applied 

simultaneously, as they are considered to he in 'this paper, they 
may also he expressed in the same form 


(E h ) 


a = k 
x c 


3;g l n 2 P 
E b 2 h 


T = k (Ea)g i nap 


where the suhcript cm signifies that the secant modulus is to he 

taken at the stress intensity = ^/cr x 2 + 3 t 2 . The values of k c 

and k: s depend on the ratio of a x to r in a manner shown in 
the appendix. 

The stress ratios are 


Bn. = 


a k ( E s) a . 

X _ c °1 

(c^x)pc ^pc ^ E s^pc 




i 




E s 


k ( E s) a . 

f_S °1 

E ps ^ E s)ps 
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For numerical calculation, a series of plates was selected, 
the first of which would buckle in the elastic range, with successive 
plates buckling at higher and higher stresses up to a stress somewhat 
beyond the yield stress. The plate material was taken as 24S—T4 
aluminum alloy with a yield stress of 46 ksi. The edges of the 
plates were assumed simply supported. Calculations for critical 
combinations of applied shear and compression were made by an energy 
method, details of which are given in the appendix. 

Results of the calculations for three plates are shown in 
figure 2. When R c and R g are used as coordinates, the curves 

of figure 2(a) are obtained. When R c and R g are modified to 

allow for variations in moduli, the curves of figure 2(b) are obtained. 
When plotted in this way the separation of the curves is appreciably 
reduced. 


The plate with the lowest value of 


jt 2 P 

b 2 h 


buckles elastically; 


the modulus ratios are therefore all unity and the corresponding 
interaction curve is the parabola of equation (1). The plate with 

2 

the intermediate value of JL - ^ buckles- elastically in pure compression 

b 2 h 

but plastically in pure shear; the corresponding interaction curve 
falls slightly below the parabola in figure 2(b). The plate with 

2 

the highest value of buckles plastically for all combinations 

_ ‘ b 2 h 

of shear and compression and its interaction curve falls still 
■farther below the parabola in figure 2(b). The sum of the abscissa 
;and the square of the ordinate (see equation (1)) is always unity 
at the ends of the interaction curves but may differ slightly in 
between. For the top (elastic) curve, the sum is always unity; 
whereas, for the bottom (plastic) curve, the sum may drop to 0.95. 

Since this difference is of the same order of magnitude as the 
inaccuracies in knowledge of material properties at very high 
stresses, the difference is not of too great consequence. As a 
practical matter it seems advantageous to set up the algebraic 

relation 


B. 


(*B> 


pc 


(E a ) 


a i 


B 


(E a ) 


ps 


s (E a ) 


= 1 


j|4ich applies to the top curve of figure 2(b). 


( 2 ) 
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By use of this relation., it is possible to discover whether 
any given combination of stresses and t will cause a plate to 

buckle. The stress intensity a j_ = \jo-y£ + 3t^ is known; therefore, 

(E q ) is known. The values of the stresses (a-y-) and T~, a 
' s v Ji, pc Ps 

under the two simple loadings can be found from references 2 and 3 ; 

at the same time the moduli (Eq) and (Eq) are also determined 

pc ps 

These quantities are substituted into equation (2) ; the plate will 
or will not buckle depending upon whether the’ value of the left-hand 
side of equation ( 2 ) is greater or less than unity. 


ILLUSTRATIVE PROBLEM 


. A simply supported plate of 24s— T4 aluminum alloy with a value 
2 

of equal to 9000 psi is considered. The stress— strain curve 

b^h -g 

of the material gives values of — as in curve A, figure 1, of 

E 

reference 2. A longitudinal compressive stress of 20 ksi is applied 
to the plate and the shear stress necessary to buckle the plate is 
required. 

From reference 2, the value of (cr x )^ c is found to be 32 ksi 

with (-^3-) = 0.9^* From reference 3, the value of t is found 

\ Jii ) ps 
' / P^ / uj \ 

to be 27.6 ksi with ( -d^-) = 0.62. - The value of R can now be 

' E 'ps 

found to be ^ = O. 625 . Equation ( 2 ) therefore gives 


0.62 

(Si 


0.625 


0.9^ 

(Si 


In order to fix the value of 


v e L’ 


a value of a ^ must be assumed. 


The value of will fall between the values of cr^ for pure 
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compression and pure shear; that is, 32 < < 48 ksi. As a first 

approximation, cr^ was taken as 40 ksi; therefore, ^ was L° un< l 

equal to 0.84, Upon substitution in the formula, R 8 is found to be 
equal to 0.74. The shear stress is" then 


T ®s T ps 

= 0.74 x 27.6 
= 20.4 ksi 

The combination of a x and T results in a stress intensity 

dj. = /cr x ^ + 3 t 2 

= ^/( 20) 2 + 3 ( 20 . 4) 2 

= 40.6 ksi 

No significant change will result if the problem is reworked with 
the new value of 0 ^ . 

CONCLUDING REMARKS 


A solution has been obtained for the buckling- stress of an 
infinitely long flat plate under combined shear and longitudinal 
compression in the plastic region. A detailed calculation for a 
number of simply supported plates made of a typical aluminum alloy 
shows that an approximate algebraic relationship, equation ( 2 ), will 
suffice for determination of buckling stress. 

Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., September 26, 1949 
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APPENDIX 


ANALYSIS 


Energy expressions in terms of Cartesian coordinates (iSy 1 ) .- 

The general expression for the net strain energy of a supported 
plate in the plastic region is given in reference 2 as formula (18). 
If only longitudinal compression o x and shear t are applied, the 

separate expressions for the strain energy and the work T 

during buckling become, in terms of the deflection w at the 
point (x’,y*) in a Cartesian system. 


V. = 


D’ 


^d 2 w \ - _ p d 2 w d 2 w , p 
L dx' 2 / 2 dx> 2 3 


d 2 w 


cbc *8y 


+ 


d 2 w d 2 w 


/ o v 

w* 


2 \ 2 

w \ 


dx ’dy 1 


dx* 2 dy'^j 
(Al) 


rn _ h 
2 


+ 2t 


dw dw 


dx'dy ' 


0X (dx ! ) dx* dy ! 

where the values of the plasticity coefficients are as follows: 


(A2) 
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jf ^ in addition., a restraining medium of rotational stiffness S is 
along one longitudinal edge of the plate , the strain energy V 2 in 

the restraint itself is taken to be., from expression ( 19 ) of 
■reference 2, 


0 r-, 


Vo = 


D*e 

2b 


U 


dy*. 


y ,= y’ 


dx 1 


(A3) 


where 


6 


_ ^Sb 


i£ y* Q is the edge coordinate. 

Ene rgy expressions in terms of oblique coordinates (x,y) .— In 

order to obtain corresponding energy expressions in oblique coordinates., 
apply the transformation 


x* = x — y sin 0 

y* = y cos 0 


> 


(A4) 


where the two coordinate systems are related as in figure 1. In the 
oblique coordinates, T, and Vg, become 


V, = Si 
1 2 


r»P 


vJO 


( C 1 - C 2 tan t ♦ 2C 3 tan 2 0 + tan^) (gf) - ^(0) 


2 ‘ 




/ C 3 it tan^lfiVc^w ¥ + ( c 3 + 2 taagA agw 
Ua 2 0 coa 2 0 A® 3 */ V 0O3 2 ^ coa 2 0 / Sx 2 9y 2 


A tan% , 4g jgg* 

\ COS p 0 CO a p 


c 2 Xd 2 * d 2 w , tan p d 2 * 
4 4- 


cos 0/dx 2 ^ 


cos 


30 dy 


2 ctocdy 


dxdy cos )6 


(A5) 
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a x + 2t tan 


+ nsfy Ie 3y 3x47 003 0 


V 2 = 


2b^ cos 0 cos 0 


+ tan 


'l 

F \d xj 


(A7) 


Energy relation at buckling .— When buckling occurs j, the energy 
relation is as follows: 


T = T 1 + y s 


Deflection surface .- In order to obtain a practical solution the 

deflection function or a reasonable approximation of it is required 
As in references 1 to 3* the deflection surface is taken to be 


W = - 1 ) + [ l + 1) COS S3, cos 

2 L 2 k \ 2 b n \ 


When this expression with its derivatives is substituted into the 
expressions for Tj and Vgj 


+ 2t tan 0^ — ( % — + i ~ — + (- - — + - cos 0 

1 k\ \120 8 jf2/ \2 n 2j 2 


p^Ci — C 2 tan 0 + 2C^ tan 2 0 + tan^ ^ + g *|-je 2 


+ /l_A_V + I +.-sk__l_ fl_i\ e 2 + /l_ji_V + l 
l 2 * 2 / 2 2b 1 3 coa V \8 U n 2p 2 


2b^X cos d 


. t (2Co + 6 tan 2 g 




€ a + (i _Jku + i 

U *2] S 
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V = 

2b23 coa 30 


^mlical combination of stre sse a » — When the foregoing values 

of Tj and V- g are substituted into equation (a8)j the following 
equation is obtained: 


a + 2t tan 0 = 


( C 1 ~ c 2 tan 0 + 2C3 tan 2 0 + tan 


Fn(0 


C3 + 3 tan 2 0 , 

+ ^ F 2 (S) 

COS^P 


Tt 2 D* 

b-L 2 h 


(■ 


— - ) COS^0 


(A 9 ) 


where 


F-,(e) = 


(jfe*S- ■!)*'*( J-a- 


F P (e) =2 


& - sti * e- - a 


1 

6 + — 

2 


+ i ~ ^y + a _ jl\ 

V-20 8 *£>; \2 


« + i 


> 


(A 10 ) 


Replacement of b^ by — 1 __ in equation (A 9 ) gives., for the 


condition of buckling. 


cos 0 
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a x + 2t tan ft 


/c^ — C2 tan j t + 2C3 tan 2 0 + tan^0 j + 


F X (€) 

/k \ 2 

l \ 


(' 


+ C 3 + 3 tan' 




0 


lt 2 P j 

b 2 h 


(All) 


When r is considered, as a given constant 3 the wave length \ 
may he adjusted to make cx x a minimum from the relation — = 0, 






which gives 


L A. i 


If- 


Fi(e) 


— C2 tan 0 + 2C 3 tan 2 0 + tan^0 


(A 12 ) 


for elastic restraint e independent of huckle wave length X„ 

Substitution of this value of in equation (All) to obtain the 

expression for a x gives 


a x = 


2 ( € ) j/c-L — Cq tan 0 + 2C^ tan 2 0 + tan^0 


f 


+ [Co + 3 tan' 


-0)f 2 ( 6 ) 


rt 2 P ! 

b 2 h 




— 2t tan ! 


(A 13 ) 


The angle } of the waves may likewise be adjusted to make a x a 


minimum from the relation 


dcr- 


< 3 (tan 0) 


= 0 } which gives 
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2 tan30 +. 2C^ tan 0 — — C 2 

C 2 tan 0 + 2C^ tan 2 0 + tan^0 


^/FjTO + 3F 2 (e)tan 0 = 


rt 2 D s 

b 2 h 


(A14) 


In the elastic range , " equations (All), \(A12), and. (Al4) reduce to 


a x + 2t tan 0 = 


(i)‘ 


COS 


+ Eikl + in ... g . |te g 0 r (t) 

^Jb'j 2 cos 2 0 


x 2 P 

b 2 h 


(x) = \/ F l (6) C °^ 2 ^ 


e) + 3F 2 (e) 


tan 0 = 


it 2 P 

h 2 h 


which agree with equations (B7), (BIO) , and (Bll), respectively, of 
reference 1. 

Computation of interaction curves.— In order to find the shear 

stress t that will cause budding for a given compressive stress cr x , 

the following procedure was used. Simply supported plates of 24s— T4 
aluminum alloy were selected for the computation; hence., 6=0, 

2-n 

Fi(e) = 1, and Fo(€) = 2. Three values of were selected; for 

bdi 

each value, equation (Al4) gives the relation between the angle 0 
and shear stress t. Pairs of values of 0 and t were then 
inserted by trial and error into equation (A13) until it was satisfied. 
Both the values of a x and the corresponding T were then divided 

by a and t^c respectively, to obtain the stress ratios. A 

iJO JJS 

complete interaction curve was obtained for each of the assumed values 
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Figure 2 shows the results of the calculation. When the stress 
ratios are used as coordinates without making allowance for changes 
in moduli., the interaction curves are as shown in figure 2(a). When 
the stress ratios are modified to allow for changes in moduli, the 
interaction curves are as shown in figure 2(h). Since the separation 
of the curves is appreciably reduced by use of the modified coordi- 
nates, curves of the type shown in figure 2(b) are recommended for 
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(a) With coordinates 
R s and R c . 


' (b) With coordinates R s and 
R c modified to allow for 
changes in- moduli. 


Figure 2.- Interaction curves for the buckling of a long flat 
plate of ZAS-TA aluminum alloy under combined longitudinal 
compression and shear. 
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